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Abstract 



The present study is based on the sample of 2.9 x 10 6 single muons observed 
by LVD at underground Gran Sasso Laboratory during 36500 live hours from 
June 1992 to February 1998. We have measured the muon intensity at slant 
depths from 3 km w.e. to 20 km w.e. Most events are high energy downward 
muons produced by meson decay in the atmosphere. The analysis of these muons 
has revealed the power index 7 of tt and K spectrum: 7 = 2.76 ± 0.05. The 
reminders are horizontal muons produced by the neutrino interactions in the 
rock surrounding LVD. The value of this flux near 90° is (6.1 ±2.7) x 10~ 13 cm" 2 
g -i sr -i The results are compared with Monte Carlo simulations and the world 
data. 
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Academy of Science, 7a, 60th October Anniversary prospect, Moscow 117312, Russia 
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1. Introduction 



The study of atmospheric muons at large depth underground is the subject of 
experimental investigations due to the following reasons. First, muons and muon- 
produced secondary particles are the background for underground detectors designed 
to search for rare events including the tasks of neutrino and gamma-ray astronomy. 
Second, the calculations of atmospheric muon and neutrino fluxes are based on the 
hypothesis about primary cosmic-ray spectrum and hadron-hadron interactions. The 
existing deep underground detectors are not able to measure muon energy for direct 
deduction of energy spectrum. But they are able to measure muon 'depth- intensity' 
curve. This curve shows the vertical muons flux as a function of the rock (water, ice) 
depth and is related to the muon propagation through the rock, muon energy spectrum 
at sea level and then, to the primary cosmic-ray spectrum. The connection to the muon 
propagation allows tests of the cross-sections of muon interactions with the rock which 
have been used in the program for the muon transport. 

The measurements of muon intensity using the single detector and technique from 
a relatively small depths to the large depth where neutrino-induced muons dominate 
are of special interest. Such an experiment observes muons at zenith angles from 
vertical to horizontal direction. The statistics for measurement of neutrino-induced 
flux in horizontal direction is small enough but the uncertainties in detecting the muon 
direction are absent. There is no need to suppress the flux of atmospheric muons by a 
factor of about 10 6 using accurate time measurements. 

The LVD structure and the complicated profile of the Gran Sasso mountains provide 
an opportunity to measure the muon 'depth- intensity' curve for slant depths from 3 
to 20 km w.e. and the neutrino-induced muon flux at horizontal direction where the 
atmospheric muon flux is suppressed due to the large slant depth. The expected number 
of horizontal events caused by neutrino-induced muons is small (about one - two events 
per one LVD tower per year), therefore for the moment we cannot make a conclusion 
about neutrino oscillations. 

In our previous paper |]J we have presented our first results on the measurement of 
the muon depth-intensity curve for the depth range of 3-20 km w.e. Since that time we 
have improved the criteria for event selection and increased the statistics. The analysis 
of Ref. was based on the events with all multiplicities. Multiple muon events 
especially for large depths are more difficult to reconstruct than single muons. To 
avoid this problem we have also performed the analysis of single muons using stronger 
criteria of run and event selection. The analysis is based on an increased statistics 
compared with our previous publications. 

In Section 2 the detector, the procedure of data analysis and conversion of the muon 
intensity to vertical are briefly described. In Section 3 the results of the analysis of the 
'depth - vertical muon intensity' relation (I^(x)) are shown. In Section 4 we present 
the analysis of neutrino-induced events. Section 5 contains our conclusions. 

2. Detector and data selection 

The LVD (Large Volume Detector) has been described extensively elsewhere |1|, 
H|, f|. The detector is located at Gran Sasso Laboratory, Italy. The minimal rock 
overburden is 3 km w.e. LVD consists of 3 towers. Each tower is made of 38 modules 
with dimensions 2.1 m x 6.2 m x 1.0 m. The data were obtained with the first LVD 
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tower from June, 1992, when it was put into operation, until February, 1998. Total live 
time was 36500 hours. The tower has the size of 13 m x 6.6 m x 12 m. Each module 
contains eight scintillator counters with active volume 1.0m x 1.5m x 1.0m and 
mass of liquid scintillator of 1.2 tons, and tracking detector attached to the bottom 
and one vertical side of the supporting structure. Each tracking detector is made of 4 
layers of tubes operating in limited streamer mode. Each layer has independent x and 
y readout strips. These established the x and y coordinates of the hits. The tracking 
system allows the measurements of particle direction with accuracy better than 0.5°. 

The mountain structure above Gran Sasso laboratory allows the measurements 
of muons which traversed slant depth from 3 to more than 12 km w.e. The depths 
correspond to the median muon energies at sea level from 1.5 to 40 TeV at zenith 
angles from 0° to 90°. In the analysis we have used the sample of events containing 
only single muons. Events with all multiplicities are usually studied in the experiments 
with cosmic rays muons. Such an analysis supposes the accurate reconstruction of 
each event. The study of depth- intensity curve with all muon events observed by 
LVD has been presented in Ref. [[J. The multiple muons have been considered as 
independent muons and the acceptance for both single and multiple muons has been 
assumed the same. This is a good and well proven approximation for the derivation 
of the all-particle primary spectrum. The tasks requires however the reconstruction of 
muon events with all multiplicities and the measurements of the direction and slant 
depths with good accuracy. This is more difficult task for multiple muon events than 
for single ones. Most muons traverse small rock thicknesses. If the slant depth for 
a small fraction of these events is defined wrong, then the intensity in this direction 
will not change much. However these erroneously reconstructed events can contribute 
significantly to the muon intensity at large depth. To be sure in the precise event 
reconstruction in this analysis we have dealt with single muons only. The size of one 
LVD tower is small enough and more than 90% of muon events are single muons. 
The number of muons in bundles is about 10-12% of the total number of muons. 
Possible uncertainties from neglecting multiple muons are less than or comparable to 
the errors from including multiple-muon events with erroneous reconstruction. In the 
case of single muon analysis we need to correct the absolute intensity for the number 
of unreconstructed events (multiple muons and muon induced cascades). 

The trigger for muon events has been defined as follows: i) energy deposition is 
grater than 30 MeV in at least two scintillator counters in two different modules; ii) 
hits in at least three layers in any three tracking detectors (hits in at least one layer 
per detector). The data runs have been selected as follows: runs have been accepted 
if they had more than one-hour duration and a counting rate is within a 15% range 
around the mean value for the set of runs. Moreover we have required at least 36 from 
38 tracking modules and 240 from 304 scintillator counters of the first tower to be 
operated in any particular run. These criteria ensured the full and uniform acceptance 
of the detector. The final muon sample after these cuts consisted of 3151580 events. 
2877659 (91%) events have been reconstructed as single muons. Multiple muons and 
muons accompanied by cascades constituted 9%. All reconstructed single muons were 
binned in two dimensional array with a cell size of 1° in azimuthal angle <fi and 0.01 in 
cos#, where 9 is the zenith angle. The accuracy of the reconstruction has been checked 
by observation of the Moon shadowing effect with single muon data []5J. It is better 
than 0.65°. 
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The acceptance for each angular bin has been calculated using the simulation of 
muons traversing LVD taking into account the detector response. The thickness of rock 
crossed by the muon was determined from the mountain map. 

The angular distribution A^(0, cos 9) obtained in the experiment has been converted 
to the 'depth-intensity relation' using the formula: 

j , x _ Y,ij N^Xrg^j, COS 9j)) 

11 m J2ij{S(x m ((j)j, cos 6i)) x e(x m (0j,cos0i)) x fiy x T) 

where the summing up has been done over all angular bins (<ftj, cos 9i) contributing to 
the depth x m ; S(x m ((f)j,cos9i)) is the cross-section of the detector in the plane per- 
pendicular to the muon track at the angle (<pj, cos#j); e(x m (4>j, cos 6^)) is the efficiency 
of muon detection and reconstruction; fi^ is the solid angle for the angular bin; T is 
the live time. Angular bins with e less than 0.03 were excluded from the analysis. The 
acceptance A is defined as follows 

A(cos 9, 0) = e(cos 9, 0) x fi(cos 9, <p) x S(cos 9, <p) (2) 

and is shown in Figure 1. 

The intensity of muons at zenith angle 9 was assumed to be related to the vertical 
intensity Iq through the relation 

I(x m ,9) = I (x mi 9 = 0°)/cos9* (3) 

where 

I?,(x m , cos9 = 1) 

is the ratio of calculated muon intensity at cos# = 1 to that at cos^j. This relation 
is valid for muons of atmospheric origin if we neglect the contribution of the prompt 
muons from charmed particles. According to the LVD data the ratio of prompt muons 
to pions does not exceed 2 x 10~ 3 at 95 % confidence level ||. 

For depth - intensity relation the bin width of 200 m. w. e. has been chosen. For 
depth more than 9 km w.e. we have chosen bins with the width of 500 km w.e. to 
increase the statistics for each bin. The conversion of muon intensity to the middle 
points of each depth bin has been done using formula: 

W=C(x m )^ (5) 

where I™(x m ) and the measured and calculated muon intensities at the 

weighted average depth x m which corresponds to the depth bin with the middle value 
of Xi\ I™{xi) and I°( the derived and calculated muon intensities at the depth 

Xi which is the middle point of the depth bin. The values of x m have been obtained 
by averaging the depths for all angular bins contributing to the given depth bin with 
a weight equal to the detected number of muons. To calculate the muon intensities 
at x m and Xi we have used the muon spectrum at sea level with previously estimated 
parameters (see also equation |^) and the simulated muon survival probabilities. 
Since the width of depth bins is quite small (200 m w.e. for depth bins with high 
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statistics) and the number of angular bins contributing to each depth bin is quite large 
(several hundreds), the conversion factor does not exceed 10%. 

3. 'Depth — vertical intensity' relation in Gran Sasso rock 

To calculate the intensity of muons underground requires the intensity of muons at 
the surface as a function of energy and zenith angle, and the survival probability as a 
function of slant depth of rock traversed: 

I,(x,cos9) = r P{E^x) dI ^ E /^ COSd) dE,, (6) 
Jo dE^o 

where P(E^ , x) is the probability of muon with an initial energy E^q at sea level 
to reach the depth x and dl^^E^, cos 9) /dE^ is the muon spectrum at sea level at 
zenith angle 9. The intensity at the surface in the units of (cm 2 s sr GeV) -1 can be 
approximated by 



dIfj,o(E f j - o, cos ( 
dE^o 



A x 0.14 x £~ 7 



1 0.054 \ 

1.1 £ M cos6>* + 1 l.lff M o cos6>* I y'l 
~r 115GeV 1 ~r 850GeV ) 

where the values of cos9 have been substituted by cos9* which have been taken from 
Ref. §. According to Ref.g cos9* = E^ K (cos9 = l)/E^ K (cos9), where E% K are 
the critical energies of pions and kaons. Equation [7| has been obtained under a simple 
assumption of scaling in the high-energy hadron-nucleus interactions. Under this as- 
sumption the power index of primary spectrum, 7, is expected to be equal to that of 
meson (pion + kaon) spectrum, 

The muons were tracked trough the rock using propagation code MUSIC [[| to 
calculate the muon survival probabilities P(E^q,x). The stochasticity of all processes 
of muon interaction with matter (nuclear interaction, pair production, bremsstrahlung 
and ionization) has been taken into account. The cross-sections were taken from Ref. 
fT0|, 0, |12||. The muon intensities calculated with bremsstrahlung cross-section from 



Ref. [12| are lower than those with bremsstrahlung cross-section from Ref. ]13[ which 



was used in our previous paper |J. Using of the cross-section from Ref. |I3] will result 



in a higher power index (softer muon spectrum) compared to the cross-section from 



Ref |]I2| . The difference in power index is of order of 0.01. 

The measured 'depth - intensity' curve is shown in Figure 2 together with the best 
fit. The underground muon flux observed at a slant depth x and zenith angle 9 has 
the two-component nature and can be presented as: 

I,(x,9) = I^(x,9) + I^(9) (8) 

where lj^'(x, 9) is the contribution of atmospheric muons and 1^(9) denotes the contri- 
bution of muons from neutrino interactions in the rock surrounding detector. For slant 
depths of 13-20 km w.e. the muons seen in LVD are of the later origin. The last ex- 
perimental point in Figure 2 corresponds to the neutrino-induced muon flux. This flux 
was measured at the depth of 13-20 km w.e. To convert the flux of neutrino-induced 
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muons to vertical intensity we used the calculated ratio of horizontal and vertical fluxes 
of neutrino- induced muons which is equal to 2.1 at energy threshold of 1 GeV for most 
models of atmospheric neutrino production: 

i;(x,e = o°) = i;(x,e = 9o°)/2.i (9) 

The 'depth - intensity' curve has been fitted with the calculated function with 
two free parameters: additional normalisation constant, A, and the power index of 
atmospheric pion and kaon spectrum, 7. As a result of the fitting procedure the 
following values of the free parameters have been obtained: A = 1.59±0.50, 7 = 2.76± 
0.05 for muon energies at sea level from 1.5 to 40 TeV. The errors of the parameters 
include both statistical and systematic uncertainties. The latter one takes into account 
possible uncertainties in the depth, rock composition, density and the uncertainty in 
the cross-sections used to simulate muon transport through the rock. These values 
are in good agreement with the results of similar analysis performed for muon events 
with all multiplicities observed by the first LVD tower during 21804 hours of live time: 
A = 1.95 ± 0.50, 7 = 2.78 ± 0.05. Note that the estimates of the parameters A and 7 
are strongly correlated. The larger the value of 7 is, the larger the normalisation factor 
A should be. 

We have repeated the fitting procedure for restricted depth ranges. The results of 
this test are presented in Table 1. The results show that the power index is the same 
within errors for all depth ranges. 

Neutrino-induced muon flux has not been included in fit procedure but has been 
added to the best fit at 2.5 x 10~ 13 cm _2 s _1 sr _1 level. Dashed curves in Figure 2 show 
possible values of muon intensities if we take into account uncertainties in the calcu- 
lation of atmospheric neutrino spectrum at sea level, structure functions, corrections 
for quasielastic scattering and energy threshold of detector. The experimental value 
(2.9±1.3) x 10~ 13 cm _2 s _1 sr _1 is in agreement with the calculated one (2.5±0.5) x 10~ 13 
cm _2 s _1 sr _1 . It also agrees with the compiled world results on underground muon in- 
tensities presented by Crouch in Ref . [|i~4H where the flux of neutrino-induced muons is 
equal to (2.17 ±0.21) x 10" 13 cm~ 2 s^sr -1 . 

If the formula from Ref. [ H| is used for the muon spectrum at sea level instead of 



eq. (0), the best fit value of 7 will be decreased by 0.04-0.05. 

The value of 7 obtained with LVD data is in reasonable agreement with the results 
of other surface and underground experiments: DEIS [TE\, MUTRON |I7| , MlPhI 



(the energies of these experiments correspond to first few points of our depth-intensity 
curve), ASD PJ, NUSEX PJ, MACRO [El, MSU [E2j (if we consider the muon spec- 



trum from |15j in the latter case). LVD data disagree with results of Baksan Scintillator 
Telescope and KGF |23| , [24| . The difference here is likely due to the different methods 
of measurements, the applied analysis procedure in each experiment and uncertainties 
in the knowledge of overburden composition. 

4. Neutrino-induced muons. 

Let us describe the evaluation of horizontal neutrino-induced muon flux in more 
details. High energy neutrinos will produce high energy muons in the rock. These 
muons will have enough energy to traverse the entire detector. Reconstructed muons 
traversed rock thickness greater than 12 km w.e. have been considered as candidates 
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for neutrino-induced muons. These depths correspond to the zenith angles more than 
85°. We have recorded 95 such candidates during 36500 hours of LVD life time. A 
careful visual scan of all these tracks eliminated 5 candidates from the sample because 
of confusion in the pattern recognition. 

Since the timing of the LVD experiment (12.5 ns) is not sufficient to determine the 
direction of a track crossing one tower there is a two-fold ambiguity in the direction for 
each measured track. In other words LVD cannot discriminate between muon direction 
(0, 4>) and (180° — 0, 180° + <fi). For < 90°it is reasonable to assume that muons come 
from above, since the rock thickness above the horizon is smaller. Gran Sasso mountain 
has a very complicated profile and for many bins at ~ 90° with x > 12 km w.e. the 
slant depth for inverse direction x\ is appear to be less than 8 km w. e. The muon 
intensity for 8 km w.e. is 80 times greater than the intensity for 12 km w.e. In this 
case we assume that the muon came from the direction with smaller slant depth. Near 
horizontal muons with reconstructed slant depth greater 12 km w.e. and slant depth in 
opposite direction less than 8 km w.e. were excluded from neutrino-induced candidates 
and considered as bins (0, 180° + 0). Totally we had 67 such events. 

Some angular bins with slant depth greater than 12 km w.e. are surrounded by 
bins with smaller slant depths. According to the calculations of Ref. || the average 
angular deviation of muons is 0.45° at 10 km w.e. and it is caused mainly by multiple 
Coulomb scattering. The probability of muon coming from the direction with smaller 
slant depth is greater. We have considered such a muon as coming from the direction 
with smaller depth assuming that it had been recorded into the bin with greater depth 
due to the reconstruction error or scattering. 

Five muons produced in neutrino interactions with surrounding rock have survived 
all cuts for slant depths x > 13 km w. e., for x > 14 km w. e. we have found four such 
events and for x > 15 km w.e. there are two neutrino-induced muons. 

A Monte Carlo has been used to estimate the expected number of neutrino-induced 
muons. The spectrum of neutrino-induced muons has been calculated following the 
formula: 

dj± = rij± x p(E^) iK (10) 

dE^ Je^ dE y dEf, 

dN v jdE v represents the neutrino spectrum, P(E U , E^/dE^ is the probability that a 
neutrino produces a muon in the interval (E^^E^ + dE^). We have used the Bartol 
neutrino flux fl2"5"f , which has a systematic uncertainty ±14%, and the Morfin and Tung 
[2~o], as well as Duke and Owens ||27| parton distributions functions, which result in 



less than 1% difference in muon spectra. The major sources of uncertainties in the 
neutrino-induced muon flux are the uncertainty in the neutrino fluxes and neutrino 
cross-sections because of required extrapolations of the structure functions to small 
x « 10 -4 . For neutrino-induced muons the calculations of atmospheric neutrino flux 
by various authors differ by as much as 17%. Different standard parameterisation of 



charged current cross sections also differ by as much as 13% |28j ). 

Table 2 shows the number of muons observed by LVD during 36500 hours as well 
as calculated values (the uncertainty of calculations is 20%). 

The observed number of muons at large slant depth agrees with predictions within 
errors. The measured neutrino induced horizontal muon flux is (6.1 ± 2.7) x 10~ 13 
cm -2 s" 1 sr" 1 , while the calculated one is (5.2 ± 1.1) x 10~ 13 cm -2 s" 1 sr _1 . Our 
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measured value agrees with the results of other experiments: Soudan-2 (5.00 ± 0.55 ± 
0.51) x 10~ 13 cm~ 2 s _1 sr _1 | 29| , Frejus (the flux recalculated for our energy threshold) 



(4.77 ± 0.86) x 10~ 13 cm" 2 s" 1 sr" 1 [||, in South Africa mine (4.59 ± 0.42) x 10~ 13 
cm" 2 s" 1 sr- 1 |3I|. 

5. Conclusions 

We have measured the underground muon intensity as a function of the slant depth 
in the range of 3 - 20 km w.e. The analysis of 'depth-intensity' relation in the depth 
range 3-12 km w.e. has been done to obtain the power index of differential energy 
spectrum of the pions and kaons in the atmosphere 7 = 2.76 ± 0.05 in the energy 
range of 1 - 40 TeV. The errors include both statistical and systematic uncertainties 
with the systematic error due to the uncertainty of the muon interaction cross-sections 
dominating. Our results are in good agreement with other experiments. Muons tra- 
versed slant depth more than 13 km w.e. were analysed to obtain the horizontal flux of 
neutrino-induced muons. This flux is equal to (6.1 ± 2.7) x 10~ 13 cm -2 s^ 1 sr _1 and is 
consistent with our calculations and reported results of other experiments. Our fit to 
this data at the 'depth-intensity' curve is in a good agreement with the fit of Crouch, 
which is a summary of various experiments. 
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Table 1: The value of power index of meson spectrum for various depth ranges. The 
errors are statistical only. 



Depth interval, km w.e. 


7 


X 2 /d.o.f. 


3-12 


2.76 ± 0.02 


25.8/34 


4-12 


2.78 ± 0.03 


19.0/29 


5-12 


2.79 ± 0.04 


16.8/24 


6-12 


2.82 ± 0.06 


14.8/19 


7-12 


2.94 ± 0.14 


11.4/14 


8-12 


2.76 ± 0.22 


6.7/9 


9-12 


2.60 ± 0.50 


3.3/4 



Table 2: The number of muons at large depths observed by LVD (N obs ) and calculated 
values (N°f m -atmospheric muons, N" - neutrino induced muons, N tot =N°f m + N°f m ). 



h, km.w.e. 


jvratm 




ATM 


N obs 


h > 13 km.w.e. 


0.45 


3.75 


4.20 


5 


h > 14 km.w.e. 


0.31 


2.92 


3.23 


4 


h > 15 km.w.e. 


0.03 


1.50 


1.53 


2 
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Figure 1: The LVD acceptance for single muons as a function of the cosine of the 
zenith angle 6 and azimuthal angle 0. The angular cell for this plot was chosen as 0.02 
(cos#) x 2° ((f)). 
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Depth In Gran Sasso rock, km w. 



Figure 2: 'Depth - vertical muon intensity' relation in Gran Sasso rock. LVD data are 
presented together with the best fit (solid curve). Dashed curves show the calculated 
intensities for maximal and minimal contributions from neutrino-induced muons (see 
text for details). 
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